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The Possibility of Kelvin-Helmholtz Instability in Solar 
Spicules 


A. Ajabshirizadeh • H. Ebadi • R. E. Vekalati • 
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Abstract Transversal oscillations of spicules axes may 
be related to the propagation of magnetohydrodynamic 
waves along them. These waves may become unstable 
and the instability can be of the Kelvin-Helmholtz type. 
We use the dispersion relation of kink mode derived 
from linearized magnetohydrodynamic equations. The 
input parameters of the derived dispersion equation, 
namely, spicules and their ambient medium densities ra¬ 
tios as well as their corresponding magnetic fields ratios, 
are considered to be within the range 0 — 1. By solving 
the dispersion equation numerically, we show that for 
higher densities and lower magnetic fields ratios within 
the range mentioned, the KHI onset in type ii spicules 
conditions is possible. This possibility decreases with 
an increase in Alfven velocity inside spicules. A rough 
criterion for appearing of Kelvin-Helmholtz instability 
is obtained. We also drive a more reliable and exact 
criterion for KHI onset of kink waves. 
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1 Introduction 


Coronal heating, the rapid increase in the tempera¬ 
ture of the solar corona up to mega kelvins is still a 
challenging problem of solar physics . Since its discov¬ 
ery in 1943 by Edlen (lEdlerJIl943 ). various theories 
have been suggested to explain this phenomenon. One 
dominant theory attributes coronal heating to the me¬ 
chanical energy of sub-photospheric motions which is 
continually transported into the corona by waves and 
dissipated somehow leading to the energy deposition 
I Zaaarashvili fc Erdelvi 2009 ). 

Acoustic waves as the prime candidate for transport¬ 
ing energy to the corona were discarded when it was 
shown that the energy flux carried by them is too 
small to have noticeable contribution to the heating 
of corona and that solar mag netic f i elds a re often 
conce ntrated in hot regions. (lAthavl Il982t iHollwea 
1990l) . So another scenario, MHD waves, which are 
supported by magnetic and density characteristics 
of solar atmo sphere were suggested as the energy 
transporters (Zaaarashvili fc Erdelvi _ 20091 Hollwea 


199C: Matsumoto k. Shibata 2010: Suzuki fc Inutsuka 


20061 ). 


Dynamic chromosphere is dominated by various fine 
features among which spicules with concentrated mag¬ 
netic fields and mass densities can be considered as 
appropriate media for MHD waves and behave as wave 
guides to carry energy from the photosphere to the 
corona. MHD waves propagating along spicules may 
become unstable and the expected instability can be of 
the Kelv in-Helmholtz type leading to arn bient plasma 
heating (|Zhelvazkov fc Zaaarashvilil 120121) . The three 
phases including energy generation by sub-photospheric 
oscillations, carrying the energy by MHD waves propa¬ 
gating along spicules and energy dissipation by Kelvin- 
Helmholtz instability are still the subjects of research. 
Spicules are chromospheric jet-like structures shoot¬ 
ing up plasma from the photosphere to the corona. 







































2 


They are observable at the limb of the Sun in chro¬ 
mospheric spectral lines mainly Ha and Can H-line. 
Since their discovery in 1877 by Secchi, physical prop¬ 


studied extensively (Beckers 

1968. I 972 I: Sterlinel 

2000 

Zaaarashvili & Erdelvil 2009 

: TsiroDOula et al. 

2012 

Camoosl 19841). De Pontieu et al. (20071) identified an- 


other type of spicules with different dynamic proper¬ 
ties and called them type ii spicules to be distinguished 
from classical or type i spicules. Although the general 
properties of type i spicules vary from one to another, 
they have diameters of about 700 — 2500 km and lengths 
of about 5000 — 9000 km. Their lifetimes lie within 5-15 
min, the temperatures are estimated as 9000— 17 000 K 
( Becker jl968 . 1972 : Zaaarashvili fc Erdelvil2009ll . The 
temperature of the ambient medium of spicules, chro¬ 
mosphere, is reported in the range of 10"^ — 10® K. 
On the other hand, type ii spicules have shorter life¬ 
times about 10 — 150 s, smaller diameters (< 200 km) 
and lengths from 1000 km to 7000 km. The short life¬ 
times of type II spicules may be either due to their 
rapid diffusion or heating and their rapid disappearance 
can be associated with the Kelvi n-Helmholtz instability 
( Murawski fc Zaaarashvili 2010ll . 

Among dynamical behaviors, oscillatory motions as 
evidences for MHD waves propagating along spicules 
have been further uncovered through imaging and spec¬ 
troscopic observations from space and ground based 
facilities. Observation of kink modes which were 
accompanied with transve rse oscillations of s picule s 
axes was first repo r ted b y lKukhiamdze,^^,a^ (12001 
Kukhianidze et al. (2006); IZaaarashvili et ~ 


by analyzing the height series oi Ha spectra in so¬ 
lar limb spicules observed their transverse oscillations 
wit h the estima t ed per i ods of 2 0 — 55 s and 75 — 110 
s. lEbadi et~^ ( 2012 1: Ebadi (2013) based on Hin- 
ode/SOT observations estimated th e oscillation period 
of spic ule axis around 90 and 180 s. IZaaarashvili et al. 
(l2007ll using doppler shift-oscillations method mea¬ 
sured the periods of spicular oscillations that could be 
due to the kink waves propagation in spicules. Es¬ 
timates of the energy flux carried by MHD (kink or 
Alfven) waves in spicules and comparisons with ad¬ 
vanced radiative magnetohydrodynamic simulations in¬ 
dicate that such waves are energetic enough possi¬ 
bly to heat the q uiet corona ( De Pontieu et al.l 12007 : 
Ebadi et al.ll2012 ). The velocities measured for spicules 
are in the range of about 20 — 25 km/s. The new type 
spicules have higher velocities of order 50 — 150 km/s. 
Hydrodynamic and magnetohydrodynamic flows are 
often subject to Kelvin-Helmholtz instability (KHl). 
In hydrodynamics KHl was first described by Kelvin 
(1871) and Helmholtz (1868). It results from velocity 


shears between two fluids moving across an interface 
with a strong enough shear to overcome th e restraining 
surface tension force (jDrazin fc Reidlll981II . In magne¬ 
tohydrodynamics the instability is generated by shear 
flows in magnetized plasmas. When the flow speed 
excee ds a critical value, KHl can arise (jChandrasekhar 
19611) . KHl in sheared flow configurations is an effi¬ 


cient mechanism to initiate mixing of fluids, transport 
of momentum and energy and development of turbu¬ 
lence. Spicules as jet-like structures flowing in ambient 
plasma environment can be subject to KHl under some 
conditions. 


ZhelvazkovI (|20121) studied the conditions under which 
MHD waves propagating along spicules become unsta¬ 
ble because of the Kelvin-Helmholtz instability. He 
considered spicules as vertical magnetic cylinders con¬ 
taining fully ionized and compressible cool plasma sur¬ 
rounded in a hot and incompressible plasma environ¬ 
ment with density 50 — 100 times less than that of 
the spicules. Considering mentioned conditions about 
spicules and their environment he derived the disper¬ 
sion equation of MHD waves which was sensitive to the 
parameters rj (density ratio of outside to that inside of 
spicule), b (the ratio of the environment magnetic field 
to that of the spicule) and Ma (the ratio of flow speed to 
Alfven velocity inside spicule) called Alfven-Mach num¬ 
ber. With r] = 0.01 and b = 0.36 and solving the disper¬ 
sion equation for sausage and kink modes, Zhelyazkov 
concluded that only kink modes can be affected by the 
KHl and the critical Alfven-Mach number needed for 
KHl onset is Mac = 12.6 corresponding to jet veloc¬ 
ity of order 882 km/s. For rj = 0.02 and b = 0.35 
the critical Alfven-Mach number was 8.9 meaning that 
jet speeds of order 712 km/s or higher are required to 
ensure that the Kelvin-Helmholtz instability occurs in 
kink waves propagating along spicules. The conclusion 
was that these speeds are too high to be registered in 
spicules. S o KHl in cylindri cal model of spicules can 
not happen) Zhelyazkov 2012 ll . 

In this paper we study the possibility of KHl in spicules 
in a wider range of input parameters 77 and b and apply 
the actual conditions of spicules. Spicules are consid¬ 
ered straight cylindrical jets flowing in ambient plasma 
environment. The basic of our investigation is the dis¬ 
persion equations of kink waves derived from linearized 
MHD equations ( Zhelvazkov 2012 : Edwin fc Robert^ 
1983ll . 


2 A summary of observational data 

Spicules physical properties as well as their ambient 
environment characteristics have been widely investi¬ 
gated. Among them densities and magnetic fields play 
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key role in KHI onset through the input parameters 77 
and b. A summary of observational data on densities 
and magnetic fields inside and outside spicules are listed 
in tables [T] and m respectively. 

It can be seen that despite providing significant data, 
there are not still definitive and certain results and even 
sometimes m ay be conflicting (|TsiroDOula et al.N2012 : 
Sterling! 200flll . To study the KHI in spicules conditions 


and according to the data in tables [T] and [2] we assume 
a wide range of values for the input parameters of the 
dispersion equation, 77 and b. This range can be consid¬ 
ered between 0 and 1 . 

Since ^ ^ we have used number density instead of 

mass density in the tables. 


as follows: 

^ + (U . V)pi + poV • vi = 0, (3) 

1 

+ Po(U • V)vi-h V(pi H-Bo-Bi) (4) 

at po 

- —(Bo-V)Bi =0, 

Mo 


^ + (U.V)pi+7MoV-vi = 0, (5) 

^ + (U.V)Bi-(Bo-V)vi+BoV-vi=0, (6) 


3 Theoretical modeling 


We consider spicule as a vertical cylinder aligned with 
the z axis of the coordinate system(z is considered per¬ 
pendicular to the surface of the Sun) with radius a filled 
with ideal compressible and fully ionized plasma of den¬ 
sity Pi and pressure Pi embedded in a plasma of density 
Pe and pressure Pe- Subindexes i and e refer to internal 
and external quantities. Inside (outside) the spicule is 
immersed in a uniform vertical magnetic field Bi {B^) 
directed along z axis. The plasma inside (outside) the 
spicule flows along the spicule axis with speed Ui {Ug). 
For convenience the ambient medium can be consid¬ 
ered as a reference frame. In this case U = Ui — Ug 
is the relative flow velocity. The equilibrium condition 
requires that the total pressure including gas and mag- 
netic pressures inside an d outside the spicule balance 
( Edwin fc Robert j 1983 1: 


Pi + Bf/2p = pg + Bl/2p, 


( 1 ) 


where p is the magnetic permeability. The density ratio 
is obtained as follows: 


^ ^ 24 + IVm 

Pi 2c2g -h -fv\^ ’ 


( 2 ) 


where Csi = \/^Pi/pi {csg = yJypTffTg) and VAi = 
Bi / y/ppi {vAg = Bg/y/ppg) are the sound and Alfven 
speeds inside (outside) the cylinder, respectively and 7 
is the ratio of the specific heats. We consider the equi¬ 
librium situation and perturb the system linearly, then: 
P = Po +Pi, P =Po+Pi, V = U-I-Vi and B = Bo-f-Bi, 
in which the indices 0 and 1 refer to the equilibrium 
and small perturbed situations, respectively. Using the 
basic equations of compressible, adiabatic, inviscid and 
ideal MHD including the continuity of mass, momen¬ 
tum, energy, induction and solenoid equations and lin¬ 
earizing them, the ideal linearized MHD equations are 


V-Bi=0. (7) 

Because the mass density of spicules in their length lim¬ 
its does not change appreciably, the gravity force term 
has been omitted in equation 01 The effects of grav¬ 
ity stratification, twist and magnetic expansion may 
be investigated in future works. By applying bound¬ 
ary conditions according to which the perturbed inter¬ 
face has to be continuous and total pressure includ¬ 
ing gas and magnetic pressures inside and outside the 
spicule must balance at the interface, and after some 
manipulations the dispersion equation of MHD waves 
for various modes propagat i ng along the spicule is ob- 


tained as (Zhelvazkov 

2012 : 

Terra-Homem et al. 

2003 

Nakariakov] 2007) 


Pi Im[Kia) 

-[(^ - k - U)^ - klv\i]ng$^^^ = 0, 


where 




and 


( 8 ) 


' Kminga) 


[(^-k.u)^-fc,^cf][(c.-k.u)^-fcM] 

(c? + - k • U)2 - klc\] 


(9) 


Ct = 


CsVA 


( 10 ) 


(cH 772)1/2- 

Im and Km are modified bessel functions of the first 
and second kind, respectively and l!^ and K'^ are their 
derivatives with respect to radial component of cylin¬ 
drical coordinates, m stands for various MHD modes 
and 777 = 1 represents the kink mode in the dispersion 
equation, n and ct are called wave attenuation coeffi- 
cient and tube velocity, respectively ( Edwin fc Robert^ 
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Table 1 A summary of mass densities and magnetic fields inside spicules 


Height(km) 

ni{cm~'^) 

BiiG) 

Source 

Investigator 

2000 - 10000 

2.2 X 10^^ - 3 X 10^^ 

18-76 

Beckers, 1968,1972 

- 

6000 

5 X lO'" - 1 X 10"^ 

- 

Beckers, 1972 


4000 - 10000 

2.11 X 10^^ - 3.5 X 10"" 

- 

Murawski & Zaqarashvili, 2013 

- 

4500 

4500 

6000 - 10000 
5000 

10“ - 10"^ 

6 X lO^'’(average) 

6 X 10^° - 1.2 X 10“ 
1.1 X 10“ - 2 X 10^“ 

6 X 10“ 

2.4 X 10“ - 4.4 X 10“ 

10-76 

10 & 30 

8-16 

12 - 15 

Tsiropula et ah, 2012 

Krat and Krat, 1971 
Alissandrakis, 1973 
Alissandrakis, 1973 

Krall et al., 1976 

Krall et al., 1976 

Kim et ah, 2008 

Lopez Ariste & Casini, 2005 
Singh and Owivedi, 2007 
Zaqarashvili et al., 2007 

Base-top 

1.3 X 10“ 

4.1 X 10“ - 5.5 X 10“ 
1.5 X 10“ 

20 

Campos, 1984 

- 

- 

3.3 X 10“ 

- 

Sterling, 2000 

- 

- 

- 

40 

Centeno et al., 2009 

Centeno et al., 2009 

- 

- 

10(quiet sun regions) 
50(active regions) 

Centeno et al., 2009 

Ramelli et ah, 2006 

- 

- 

18 - 76 & 12 - 56 

Yeon-Han Kim et ah, 2008 

Yeon-Han Kim et ah, 2008 


Table 2 A summary of mass densities and magnetic fields outside spicules (chromosphere) 


Region 

ne(cm 

Be{G) 

Source 

Investigator 

cool corona 

1 X 10^ 

10 

Aschwanden, 2005 

- 

2000 - 10000 km 

8 X 10^ - 5 X 10® 

- 

Aschwanden, 2005 

Fontenla et ah, 1990 





Abriel, 1976 

- 

10^ 

- 

Zaqarashvili & Erdelyi, 2009 

- 

quiet sun region 

1 X 10® - 2 X 10® 

- 



coronal holes 

5 X 10^ - 1 X 10® 

- 

Aschwanden, 2005 

- 

coronal streamers 

3 X 10® - 5 X 10® 

- 



active regions 

2 X 10® - 2 X 10® 

- 



- 

10“ - 10^® 

- 

Aschwanden et al., 2002 

- 



























5 


1983tl . Normalizing all velocities in dispersion equation 
[8] to VAi yields more appropriate form of the dispersion 
equation. 


- [{V,n - Ma? - l]«ea (11) 
KmiHeO) 


in which Vph = VphjvAi = ^/(kzVAi) is dimensionless 
phase velocity, K = k^a is wave number normalized to 
spicule radius, Ma = U/vAi is Alfven-Mach number, 
T] = Pel Pi is the density contrast and b = Be/Bi is the 
ratio of the exterior magnetic field to the interior. Con¬ 
sidering spicule as a cool plasma jet and the environ¬ 
ment as a hot and incompressible plasma leads conse¬ 
quently to Cs 0 and Cs —>■ oo, respectively and corre¬ 
spondingly the internal and external wave attenuation 
coefficients change into Kia = [1 — {Vph — 
and /Cefl = K. Using these approximations the disper¬ 
sion equation El can be more simplified and for kink 
mode (m = 1) it takes the form 


(lyt - i.=)ii - (Vi. - 


-[{Vph - Ma^ - 1] 


K[{K) 

Ki{K) 


= 0 . 


( 12 ) 


By solving equation [TT] we can obtain phase velocity 
as a function of wave number Vph{K) in which the real 
part Re{Vph) represents dispersion curve of kink waves 
and the imaginary part Im{Vph) represents the growth 
rate of unstable kink waves propagating along spicules. 


4 Numerical results and discussion 


First we sol ved equation [TT] using Muller method 
(jMnlleii Il95a l at various values of the input parame¬ 
ters r] and b. Based on the observational results, they 
both are considered to be within the range 0 — 1. In 
each case varying Alfven-Mach number Ma from zero 
to some reasonable numbers, yields instability threshold 
value called critical Alfven-Mach number Mac- When 
Ma reaches the critical value Mac, the system becomes 
unstable and the instability is of Kelvin-Helmholtz 
type. The beginning of KHI is accompanied by ap¬ 
pearing growth rates curves in the imaginary part of 
Vph{K) and changing the closed dispersion curves into 
nearly smooth ones. For input parameters b = 0.25 and 
T] = 0.1 we obtained Mac = 3.93. The dispersion (up¬ 
per panel) and instability growth rates (lower panel) 
curves for these parameters are depicted in Figure [TJ 



k.a 



Fig. 1 Dispersion curves of kink waves propagating along 
spicules (upper panel) and growth rates curves of unstable 
kink waves (lower panel) at the input parameters b — 0.25 
and rj = 0.1. 
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The results of solving equation [TT] at various values 
of r] and b within the range 0—1 are shown in Figure 
[2] The values of MacS are demonstrated by colors. 
We note that smaller values of MacS are concentrated 
in the regions with lower density contrasts and smaller 
magnetic fields ratios. These regions are specified by 
dark blue color. On the other hand red-colored areas 
present larger values of Mac- These areas correspond 
to the smaller r^’s and larger b’s. 

Figure [3] and Figure 0] show the behavior of Mac 
with respect to the parameters r] and 6, respectively. 

According to these figures it can be clearly empha¬ 
sized that Mac becomes smaller while r] increases and 
b decreases. To discuss the subject more physically, we 
employ the relation Ma = ■ It is obvious that when 

the system is in critical situation, this relation can be 
written in the form Uc = VazMac where Uc is critical 
flow speed, the minimum speed at which the flow be¬ 
comes unstable. It can be concluded that at smaller 
values of Mac, lower critical flow speeds are needed for 
KHI onset. Since the range of type ii spicules velocities 
are estimated to be between 50 km/s and 150 km/s, 
we are interested in the regions oi b — rj plane (we call 
them KHI regions) where MacS take the values that 
their corresponding C/c’s can overlap spicules speeds. 
On the other hand we know from the aforementioned 
relation that Uc is proportional to VAi too, therefore 
to determine critical flow speed of a spicule at a given 
point {b,r]), it is necessary to specify Alfveen velocity 
Vaz inside it. Mathematically Vaz is a function of den¬ 
sity and magnetic field, therefore one can find a lot of 
choices for it at the point {b,r]). To avoid the com¬ 
plexity due to the multiple choices for Vaz, we prefer to 
continue our discussion using various values of Alfven 
velocities. Considering Vaz = 40 km/s, the minimum 
and maximum values of critical flow speed can be ob¬ 
tained as Ucmin = 60 km/s and Ucmax = 625 km/s cor¬ 
responding to 77 = 1, 6 = 0.05 and 77 = 0.01, b = 1, re¬ 
spectively. The comparison between spicules observed 
velocities, (50—150) km/s, and their critical flow speeds 
,(60 — 625) km/s shows a partly overlapping. The over¬ 
lapping region (KHI region) for Vaz = AOkm/s is the 
whole area above dashed line A in Figure 0] For the 
parameters b and 77 which their values are out of KHI 
region, the critical flow speeds are higher than those 
ever observed in spicules and so can be discarded. As 
a result we can deduce that spicules which their in¬ 
put parameters (6,77) are inside the KHI region and 
flowing with the minimum speed (Uc) assigned to that 
point according to the considered Alfven velocity, can 
be subject to Kelvin-Helmholtz instability. By increas¬ 
ing Vai overlapping between spicules observed and crit¬ 
ical speeds gradually decreases and the KHI region be- 



b-il/ik 


Fig. 2 Minimum Alfven-Mach numbers in which spicules 
as cylindrical plasma jets can become unstable of Kelvin- 
Helmholtz type based on various values of 77 and b. The 
colors identify Mac’s. The determined regions correspond 
to the range of type ll spicules speeds in various Alfven 
velocities. 


M 


— Mil 



0 Oil 02 M 04 06 06 07 06 06 i 


n*p./Pi 

Fig. 3 The behavior of Alfven-Mach number in terms of 
the input parameter 77 . 
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Fig. 4 The behavior of Alfven-Mach number in terms of 
the input parameter b. 
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comes smaller indicating that the possibility of occur¬ 
ring KHI diminishes. In Figure[2]KHI regions for Alfven 
velocities of 50 km/s, 70 km/s, 80 km/s, 90 km/s and 
100 km/s are indicated with areas above dashed lines 
B, C, D, E and F. At VAi = 101 km/s and higher the 
minimum critical flow speed exceeds the maximum type 
II spicules velocity (150 km/s) leading to zero overlap. 
In other words for Vai > 100 km/s the possibility of 
KHI in spicules vanishes. 

A rough criterion for the occurrence of KH instability 
of kink waves was suggest ed by Andries & Goossens 
(|Andries fc Goossensll2001ll . In our notation it is in the 
form of 


-^AoU-b/yTJ- (13) 

We derived another rough criterion for KHI onset of 
kink waves. Since Bessel functions and their deriva¬ 
tives in the general and simplihed equations ITT] and [T^ 
respectively are appeared as the ratios, we realized that 
these ratios are real valued and moreover 

I'i{Kia) _ Ke K[{K^a) 

hinia) K, Ki{Kea)' 

By these assumptions the resulted criterion is: 

>y^(62 + l) 

which seems to give more reliable and exact predictions. 


(14) 

(15) 


KHI. According to the relation Uc = Vai Mac and sub¬ 
stituting VAi = 40 kmjs we can determine a region in 
the b — rj plane where the critical flow speeds overlap 
with the type ii spicules speeds. Out of this region UcS 
are higher than the observed speeds of spicules and 
then they can be discarded. Gonsidering larger values 
for Vai makes KHI region smaller indicating that the 
possibility of KHI onset in type ii spicules conditions 
becomes smaller. For Alfven velocities of 89 km/s and 
lower the onset of Kelvin-Helmholtz instability in type 
II spicules conditions is possible. This possibility de¬ 
pends on the area of the KHI region. The curves of 
Mac in terms of rj and b show that the threshold value 
of Alfven-Mach number becomes smaller while increas¬ 
ing r] and decreasing h. 

Solving the equation [TT] in which spicules and their 
surrounding environment are considered in more re¬ 
alistic conditions and the exerted approximations are 
removed, led to critical Alfven-Mach numbers smaller 
than those resulted from equation [ 12 ] showing that ap¬ 
pearing of KHI gets more possible. 

To provide more reliable predictions for KHI onset of 
kink waves, we derived a criterion using Bessel func¬ 
tions properties. It takes the form 

Mac > ^^( 6 ^ + 1 ). 

Acknowledgements This work has been supported 
by the Research Institute for Astronomy and Astro¬ 
physics of Maragha (RIAAM), Maragha, Iran. 


5 Conclusion 

We considered spicules as cylindrical jets shooting up 
plasma from the Sun’s surface along the z axis of the 
cylinder. Both spicules and their environment are em¬ 
bedded in magnetic fields directed in z axis. Using 
linearized MHD equations we obtained the dispersion 
equation of kink waves propagating along the spicules. 
The key parameters of the dispersion equation which 
affect Kelvin-Hehnholtz instability occurrence, are 77 , b 
and Ma- They are dehned as rj = PejPi, b = Bg/Bi 
and Ma = U/VAi, where rj and b are mass densities and 
magnetic helds ratios of outside to inside spicules and 
Ma is the ratio of spicule streaming speed to Alfven ve¬ 
locity inside it. Observational values of physical prop¬ 
erties pe, Pi, Be and Bi reported in the literature are 
employed to determine the key parameters p and b. It 
was shown that they can range from 0 to 1. We solved 
the dispersion equation at various values of 77 and b 
within the range mentioned, and obtained a critical 
Alfven-Mach number (Mac) for each ( 6 , 77 ) pair. The 
final result is an array of MacS representing the thresh¬ 
old values of Ma's in which spicules can be subject to 
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